Quantum Control Theory for State Transformations: 
Dark States and their Enlightenment 
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For many quantum information protocols sucii as state transfer, entanglement transfer and en- 
tanglement generation, standard notions of controllability for quantum systems are too strong. We 
introduce the weaker notion of accessible pairs, and prove an upper bound on the achievable fidelity 
of a transformation between a pair of states based on the symmetries of the system. A large class of 
spin networks is presented for which this bound can be saturated. In this context, we show how the 
inaccessible dark states for a given excitation-preserving evolution can be calculated, and illustrate 
how some of these can be accessed using extra catalytic excitations. This emphasizes that it is not 
sufficient for analyses of state transfer in spin networks to restrict to the single excitation subspace. 
One class of symmetries in these spin networks is exactly characterized in terms of the underlying 
graph properties. 



I. INTRODUCTION 

The holy grail of quantum information is the realiza- 
tion of a universal quantum computer, capable of solving 
a variety of tasks, including some that cannot be solved 
efficiently by classical computers. Although progress has 
been made towards this goal, it remains a daunting task. 
A more practical alternative is to develop utility modules 
capable of specific tasks such as quantum wires for state 
transfer [l|, Q or entanglement generation 043 j which 
are easier to design and manufacture. These tasks can be 
achieved with minimal or no control. Perfect state trans- 
fer, for example, can be achieved by propagation along a 
spin chain with a suitably engineered Hamiltonian [l|, yj . 
However, these schemes have obvious drawbacks such as 
the requirement for precise Hamiltonian engineering, and 
the fact that they are limited to performing a single fixed 
task. We would like to relax the constraints slightly and 
see how much is gained by adding a small amount of 
control. Do we enlarge the range of tasks that can be 
achieved and the systems that can be used? 

Recent work suggests that it is indeed possible to vastly 
increase the potential of many systems with very lim- 
ited control. Control of a single coupling and some lo- 
cal control of the end spins, for example, can be suffi- 
cient to efficiently realize a universal quantum computer 
@. The problem of realizing quantum information 
primitives such as state transfer and entanglement gen- 
eration using a minimal degree of control, and the re- 
lated task of identifying the key properties of a system 
such that those control techniques can then be applied, 
have also been studied, in particular for spin systems. 
Most studies have focused on special geometries, partic- 
ularly chains and quasi-one-dimensional systems [ll-[Ei0~ 
[§], dual-chains [10] and rings [11.], but there have been 
some studies more general networks [H, [11]. Each of 
these investigations has proved a certain information pro- 
cessing property, but a general framework for determin- 
ing the capability of an arbitrary network is still lacking. 

In this paper, we propose a unifying framework for de- 



scribing certain state transformation tasks that involve 
simultaneously transforming a fixed set of initial states 
into a fixed set of output states. This allows us to simul- 
taneously discuss processing tasks such as perfect state 
transfer and entanglement generation, rather than tackle 
each case separately. These tasks require far less than full 
control of the system, but the ability to perform them is 
restricted by symmetries of the system. In Sec.|ll]we dis- 
cuss the notions of full controllability, and introduce the 
much weaker notion of sets of accessible pairs, and de- 
scribe the role of Lie algebra symmetries, deriving bounds 
on the fidelity with which a given state transformation 
task can be realized based on symmetry properties of the 
system and control Hamiltonians. In Sec. IIIII wc demon- 
strate how the upper bounds on a state transformation 
task can be achieved within a large class of spin net- 
works by controlling a single coupling between two spins. 
This includes an exact description of one class of symme- 
tries for these systems, a necessary and sufficient condi- 
tion for their presence being that the underlying network 
topology should correspond to a bipartite (two-colorable) 
graph. In Sec. II VI it is shown how symmetry restrictions 
can be overcome by catalytic excitations, which empha- 
sizes a marked difference in state transformation capabil- 
ities in comparison to the standard assumptions made in 
the analysis of state transfer in networks [Sli] . Finally, in 
Sec. |Vl we discuss how many of the important parame- 
ters for the transfer can be estimated from experimental 
data if they are not previously known. 



II. SYMMETRIES & ACCESSIBLE STATES 

A. Controllability 

Among the most elementary concepts in control theory 
are reachability and controllability. Given a controlled 
dynamical system, a state xi is reachable from an initial 
state Xo if there exists a control sequence and some time 
T such that the corresponding trajectory of the system 
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satisfies x(0) = Xq and x(T) = Xi, and the system is con- 
trollable if any state in the state space is reachable from 
any other state Apphed to bihnear control systems 
evolving on a Lie group, one can derive explicit criteria 
for controllability in terms of the Lie algebra formed by 
the dynamical generators. In particular, this is the case 
for decoherence-free quantum systems governed by the 
Schrodinger equation il{t) — —iH{t)U{t), where U{t) is 
a unitary operator on the Hilbert space H of the system, 
and H{t) is a Hamiltonian that depends linearly on a set 
of (real- valued) controls fm{t) 



H{t) = Ho 



fm{t). 



(1) 



Here Hq is the system Hamiltonian and the Hm 
are Hamiltonians governing the interaction with time- 
varying control fields fm(t)- iHo and iHm are the gener- 
ators of the dynamics and they generate the dynamical 
Lie algebra L, which consists of all linear combinations 
and iterated commutators of these generators. For a sys- 
tem with Hilbert space dimension n the generators iHm, 
m = 0, 1, . . ., are anti-Hermitian matrices and thus can, 
at most, generate a subalgebra of the Lie algebra u{n) of 
anti-Hermitian nxn matrices. A necessary and sufficient 
condition for the system ([T|) to be fully controllable then 
is the Lie algebra rank condition L ~ u{n). If it is 
satisfied, then it can be shown that for any sufficiently 
large T, there exists a set of controls fm{t) such that any 
desired unitary gate t/targ G can be realized, i.e., 

given any ?7targ G ii{N), we have jl5j | 



C^targ ^U{T) = cxp+ < -i 



[ Ho+J2fm{t)Hmdt\ , 



(2) 

where exp_,_ denotes a positive time-ordered exponential. 
The reachable set K = exp(L), known as the dynami- 
cal Lie group is, in this case, the entire unitary group 
Applied to an iV-qubit system with a control- linear 
Hamiltonian ([IJ , this implies that the system is fully con- 
trollable if and only if L = u(2^), and the corresponding 
dynamical Lie group K = il(2^). In practice, this con- 
dition can be slightly relaxed as we can usually neglect 
the global phase of the system, in which case L = su{n) 
suffices, where su{n) is the Lie algebra of trace-zero anti- 
Hermitian nxn matrices. In this paper we consider sys- 
tems which are not fully controllable, for which L need 
not be u{n) or su(n). 



given pair is not accessible then we quantify how well the 
system can achieve the desired transformation by calcu- 
lating the fidelity of the optimal output state IV'out), 

^i-K^outlV-targ)!'. 

Accessible pairs have F ~ 1, whereas a state is dark 
with respect to a given input if _F — 0. One can ex- 
tend this concept to sets of accessible pairs. We say 
a system has a set of (simultaneously) accessible pairs 
{(K), IV'Lg)), (IV'fn), iCrg))} if the state |0 can be 
transformed into iV'taig) some set of manipulations, 
and the same set of manipulations, in the same time, 
also transforms iV'fn) ii^to IV'targ)- 

This formalism provides a simple way for describing 
state transfer. Consider two nodes of a network A and 
B. State transfer from A to B can be achieved if 

{(10)^10), |^)|0)s),(|1)a|^),|^)|1)b)} 

are accessible pairs for any arbitrary states |0) and \ip) 
as it implies 

ia\0)A + PIDaM ^ W{a\0)B + 

for any superposition a|0) -I- (In this paper, we will 
only consider the decoherence-free case, which permits 
us to use linearity of the transformation task). Previous 
studies have focused almost exclusively on systems with 
excitation-preserving Hamiltonians because this vastly 
simplifies the task. As jO)**^ is an eigenstate of such 
a system, setting \(p) — — jO)'*^""'^ shows that we 
simply require that the system has an accessible pair 
{\1)a\0)^^~\ |0)®^-^|1)b). We also mention that if we 
have an initial state -^(|00) + |11))^.^'|(/)), where A' is an 
ancilla qubit not part of the system, then the same acces- 
sible set as for state transfer allows the entanglement to 
be distributed between A' and B, ^(|00) -I- |11))s,A'|^/')- 
The concept of accessible pairs is a very weak property 
compared to full controllability. It requires only that we 
can implement one unitary operator that simultaneously 
maps the input states |V'in) to the target states |i/'^^arg)' 
as opposed to an arbitrary unitary operator. In the ex- 
ample above, this means we only need to generate one 
unitary operator U that maps |10 • • • 0) to |00 • • • 01) to 
swap the quantum states of the first and last qubit. A 
one-dimensional Lie algebra may suffice for these pur- 
poses Jj]. 



B. Sets of Accessible Pairs 

The main task in which we are interested throughout 
this paper is whether a system described by Eqn. ^ 
can be caused to transform its initial state iV'in) hi^to a 
given target state I^Atarg) through suitable manipulations 
of the control fields fmit). If this is the case, we say 
(IV'in), IV'targ)) IS an acccssiblc pair of the system. If a 



C. Reducibility and commuting symmetries 

The most important reason for lack of full controlla- 
bility of a control system on a Lie group are symmetries 
restricting the dynamical Lie algebra and dynamical Lie 
group. The relation between symmetries in the Hamilto- 
nian of the system, irreducibility and loss of full control 
has been studied in recent papers by Polack et al. 
and Sander et al. p7| . Lie algebra symmetries can be 
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classified in various ways. We shall distinguish between 
symmetries that imply reducibility of the dynamical Lie 
algebra and other symmetries. 

A dynamic Lie algebra L is said to be reducible if it can 
be written as the direct sum of simple Lie algebras, e.g., 
A necessary and sufficient condition for re- 
ducibility of L is the existence of a non-trivial Hermitian 
symmetry operator J that commutes with every element 
of L, 

[iJ„„J]=0, Vm = 0,l,2,.... 

We shall refer to such J as commuting symmetry opera- 
tors (CSOs). The commutation relations imply that the 
Hamiltonians leave non-trivial subspaces of the Hilbert 
space Ji invariant, and thus there exists a basis with re- 
spect to which J and all iJ„j have the same non-trivial 
block-diagonal structure, and we can write the Hilbert 
space ?^ as a direct sum of mutually invariant subspaces 
Ti = (B^^iHd, where K is the number of independent 
blocks. Every Hm has an eigenvector decomposition such 
that every eigenvector can be written as a linear super- 
position of the eigenvectors of J drawn from a given de- 
generate eigenspace. The evolution on each subspace is 
described by its own Hamiltonian and thus the projection 
of a state onto each subspace, ||ndV'(i)||7 is a constant of 
motion. Therefore, a necessary condition for a system to 
have an accessible pair (lipin), IV'targ)) is that the initial 
state \ipin) and target state IV'targ) have the same projec- 
tions onto each invariant subspace. 

An important subclass of CSOs are the permutation 
symmetries, corresponding to the interchange of indis- 
tinguishable physical qubits of the A^-qubit system. This 
particular subclass was referred to as an outer symmetry 
in [l7|. A general procedure to find CSO is described in 
Appendix 



D. Irreducibility and other symmetries 

Irreducibility of the Lie algebra on an n-dimensional 
subspace Hk alone, however, is insufficient to ensure full 
controllability on this subspace. Assuming the Hamilto- 
nians Hm are traceless, which can be achieved by sub- 
tracting a multiple of the identity if necessary, all irre- 
ducible components of the Lie algebra must be simple 
subalgebras of su(n), which restricts to candidates of the 
form sp{k), so{k), su(fc) (for some A: < n) or one of the 
exceptional Lie algebras. In this case, although we can 
always produce a unitary evolution that takes an initial 
state to a final state with a non-zero overlap with our 
target state, we may not be able to produce arbitrarily 
high overlap with the target [l^. For a subspace of di- 
mension n this occurs precisely when the dynamical Lie 
algebra is neither su(n) nor sp(^) 18 1. 

Example 1: A simple example is a system with 



Ho 




Hi 




The Lie algebra generated by {iHq, iHi} is a unitary rep- 
resentation of so(3) and indecomposable, but the system 
possesses an orthogonal symmetry; both matrices anti- 
commute with 




Let us define the basis to be |±A) and where |±A) are 
the eigenvectors of Hq with eigenvalues ±1 respectively. 
Hi couples the state |1) equally to the states |±A), which 
are not directly coupled. Imagine we start from the state 
|1) and wish to produce a state The symmetry pre- 
vents us from perfectly producing |^) unless 

|(A|*)| = |(-A|*)|. 

To see this, consider the unitary basis change 

1 




which maps the generators iHo and iHi onto real skew- 
symmetric matrices. Hence the evolution under any con- 
trol sequence 



B exp I < — i 



^ ^Ho + fit)Hidt^B^ 



must be described by a unitary matrix which is real. 
Since B\l) is a real vector (up to a global phase), the 
output -B|^') must be a real vector {a,b,c)^. Hence, in- 
verting B shows that we have (iAI*) = ^{a + cTV^ib) 
which, because each of the coefficients is real, have equal 
magnitudes of amplitude. 

Thus, given a maximal decomposition of the Hilbert 
space and Lie algebra into irreducible subalgebras, the 
next step is to identify symmetries such as orthogonal or 
symplectic symmetries, which are characterized by the 
existence of an orthogonal or symplectic symmetry oper- 
ator J such that every element x in the irreducible Lie 
algebra Li satisfies 



Jx = 0. 



(3) 



If all the dynamic generators of the system are of the form 
X = iHm, where Hm are real-symmetric matrices, then 
(iHm)^ J+J{iHm) = simplifies to the anticommutation 
relation 



{Hm,J}=0 Vm==0,l,2, 



(4) 



and we shall therefore refer to such J as anticommut- 
ing symmetry operators (ASOs). Given the Hamiltoni- 
ans Hm, such symmetries can easily be determined, as 
described in Appendix [X] 
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1. Simultaneous Symmetries 



sis, we have 



A given system may have several CSOs that do not 
commute with each other, meaning that block-diagonal 
structures are not simultaneously realizable. We will al- 
ways choose to analyze the symmetries that are most 
natural for the problem at hand, but one might worry 
that other non-commuting symmetries, be they CSOs or 
ASOs, could come into play. We will briefly justify that 
they do not, and therefore once we have picked a set of 
CSOs, we only have to find other CSOs or ASOs that act 
on the reduced space. 

Imagine that we have found a CSO such that the sys- 
tem Hamiltonian and the control Hamiltonians all de- 
compose into the same subspace structure, 

and that our initial state only has support on T-Li. 
Clearly, the evolution of that initial state is only de- 
termined by Hamiltonians restricted to Hi; it does not 
matter what the Hamiltonians are on the complement to 
that space. So, let's say we have now found an M which 
could either be a CSO or ASO, but does not respect the 
subspace structure i.e. it couples between the two. The 
existence of this M must depend heavily on the structure 
of the Hamiltonians on the complement space - if they 
were different, M would not exist, and then it certainly 
couldn't restrict the dynamics within "Hi. Since we've 
argued, however, that our evolution of a state supported 
on "Hi should not be affected by changing the Hamilto- 
nians on the complement space, it must be that M does 
not affect the dynamics. 



E. Maximum information transfer fidelity 

Combining the knowledge of how the system Hamilto- 
nian decomposes under basis change with the knowledge 
of the symmetries satisfied by each subsystem allows us 
to characterize what properties of the system's popula- 
tion are conserved under controlled evolution. This infor- 
mation can be used to calculate the maximum achievable 
overlap of an initial state with a target state as a measure 
of a network's suitability for certain information trans- 
mission and state preparation tasks. 

First, let {|0in)} be an orthonormal set of input states 
(we will typically take these to be determined by the 
subspace structure induced by the CSOs), which we can 
use to decompose the input state 



l^in) 



Now suppose we can implement a particular unitary op- 
erator U that maps this set to a set of output states 



|V'targ)=E/3n|Cu{) 

n 

whereas U can only transform to 

lCt) = E"»l'^oul), 

n 

giving a maximum overlap with the target state of 



(5) 



If i^' = exp(L) is the dynamical Lie group of the system, 
L being its dynamical Lie algebra, then the maximum 
transfer fidelity is 



F = maxueK KV'^tlV'targ) 



(6) 



The latter expression is not easy to evaluate in general, 
but we can give upper bounds based on our knowledge 
of the symmetries of the system. 

For example, suppose we know that the system decom- 
poses due to the presence of CSOs, i.e., H = 0^ Hd with 
dim Jid — Nd and = N, and the dynamics on each 

subspace Hd are independent. We have 

|V'in)=End|^in)=E"'il'^in^) 
d d 

l^targ) = ^n^lVtarg) = ^ ^^^^ • "^oul ) 



where ad = ||H<i?/'in|| , l^d = Hn^Vtargl! and = 

Q;;^^nd|V'in), Ifl^iut) = ^n^lV'targ) are the respective 
(normalized) projections onto the subspaces. As any vec- 
tor in the subspace Hd must remain in this subspace and 
the norm is preserved under unitary evolution, the best 
we can hope for is to simultaneously map the normalized 
for all d, in which case we obtain 



states l^ln'') o'^to 



F=\Ed<M 



(7) 



Expressing the target state in this ba- 



which is equivalent to This bound is guaranteed 

to be attainable only if all the relative phases between 
subspaces can be controlled, i.e., if the dynamical Lie al- 
gebras on the subspaces are u{Nd) or sp{Nd/2) © u(l). 
In practice, if the subspace Hamiltonians have zero trace 
then we do not have (subspace) phase controllability, and 
the most we can actually hope for is to map each 
to a state e"^'*|^Q^-{.), in which case we obtain the modi- 
fied bound F = \ J2d<^dPde''f'''\ and the original bound 
may be unattainable. Similarly, this may be the case if 
the subspace dynamics is constrained by further ASOs. 
Nonetheless the bound is often attainable for certain ini- 
tial and target states even in spite of these symmetries 
and phase constraints. 
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III. APPLICATION TO SPIN NETWORKS 

To demonstrate some of these concepts, we take the 
well-known class of examples as defined by the XX- 
coupled iV-spin networks. 

Define a graph Gq = (Vo,i?o) of the spin 
network, where Vb is the vertex set (of N ele- 
ments) of physical qubits and the edge set Eq = 
{(i,j)|qubit i interacts directly with qubit j}. The sys- 
tem Hamiltonian for the spin network is thus 



^ df^{X,X,+Y,Y,), 



(8) 



where X, Y and Z are the standard Pauli matrices and 
dij are real-valued coupling constants. The Hamiltonians 
are therefore given by real-symmetric matrices. 

We can similarly add control to the network by modi- 
fying a subset of the couplings defined by the edge set of 
a second graph, Gc — (Vb,-Ec) 



{h3}&Ec 



{X,X,+Y,Y,) 



(9) 



The controlled coupling can be manipulated either by a 
simple on-off sequence or more complicated pulses. It 
can easily be shown that the system Hamiltonian com- 
mutes with the total spin operator: [Hq, Sz] = for 
Sz '■= J2iLi s-nd thus the number of excitations in 
the system is preserved, i.e., the dynamic Lie algebra is a 
direct sum of N subalgebras, each corresponding to the 
evolution of fc G {1, . . . , N} excitations along the chain. 
The control He also preserves the coupling isotropy in 
the a:y-plane and thus also commutes with the total spin 
operator, [Hc,Sz] = 0, preserving the decomposition 
into excitation subspaces. The system is therefore not 
controllable in a conventional sense and would be unsuit- 
able for A^-qubit quantum computation without the ad- 
dition of some control that does not commute with Sz in 
order to break this symmetry and couple the sub-blocks 
of the Hamiltonian. However, for information transfer 
and state manipulations, excitation preservation is a pos- 
itive property to exploit, automatically confining the dis- 
persion of the excitation(s) within a lower-dimensional 
subspace. This will allow us to show how a single excita- 
tion can be forced to propagate to a chosen target state 
of the network. 

Within the first excitation subspace, we use the ba- 
sis \n) to denote an excitation localized on the n' 
qubit. Similarly, in the second excitation subspace we 
use 1 71, m) to denote a pair of excitations on qubits n 
and TO. In a slight abuse of this notation, given a 
state It/i) = J2n=2 '^nl'n) , we will write to denote 

J2n=2 an\l,n). 

We shall use the matrices A and C to denote the re- 
striction of Hq and Hq to the single excitation subspace. 
A is an N X N matrix with elements c?|°'' ^ for any pair 



th 



(iyj) which constitutes an edge of the underlying graph 
(and similarly for C using dfj instead). For simplicity, we 

shall take > 0. If all d\^^ — 1, then A corresponds 
to the adjacency matrix of the graph Go. Note that all 
diagonal elements of both matrices are zero. 

Although the following analysis could be conducted 
for more general controlling interactions, in the spirit of 
analyzing the capabilities of systems with minimalistic 
control, we consider controlling interactions of pendant 
type, i.e. in graph Go, vertex 1 is unconnected (meaning 
Aij = Aji =0 for all j), whereas the only edge in Gc is 
between vertices 1 and 2. Without loss of generality, we 
can take di2 = 1, corresponding to a control Hamiltonian 



He ■.= ^{XiX2 + YiY2). 



(10) 



In the first instance, we assume that the coupling 
strengths are all known. 

We first prepare the network in the zero-excitation 
state using a suitably strong background magnetic 
field 33|, and initialize a state transfer or a single ex- 
citation state preparation task by introducing a single 
excitation on the pendant qubit. Given the initial state 
|1) on the pendant qubit and a target state \tp) in the 
single excitation subspace, what is the optimal fidelity 
that can be achieved through modulation of the control 
G? 

We start by applying the symmetry conditions to the 
particular case of a pendant controlled spin network, re- 
stricted to the first excitation subspace. Let M be a 
symmetry operator of the system. Then 

MA ± AM = 0, 
MC±GM = 0, 

where the sign determines whether the symmetry is of 
ASO or CSO type and must take the same value for both 
equations. 

If M is a CSO, then |1) must be an eigenstate of M. 
Also, the eigenvectors of G, |1) ± |2) must be eigenstates 
of M. Combining these statements, |2) is an eigenstate 
of M, and this is in the same eigenspace as Any 
eigenstate of A that is in a different eigenspace of M 
must satisfy (1|G|A) = 0, and cannot be accessed. This 
allows us to show that degeneracies of the system are 
necessarily associated with symmetries. If A has degen- 
erate eigenvalues then we have an eigendecomposition 
A = J2k^kJlk, where H^ are projectors onto the re- 
spective eigenspaces, and we have either Hfe|2) = 0, in 
which case there is no overlap of state |2) = G|l) with 
the fcth eigenspace, or we can define an eigenbasis such 
that \Xk,i) = Hfe|2). As the projection of a vector onto a 
subspace cannot have more than rank 1, all other eigen- 
vectors \Xk/) in the kth eigenspace must have zero over- 
lap with state |2), and thus each \Xk,i){Xk,i\ for I > 1 
commutes with both A and G and is therefore a CSO. 
This means that when we work under the assumption of 
having found all CSOs, we can take the subspaces to be 
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non-degenerate. Henceforth, we shall use T-La to denote 
the smallest subspace with support on |2), i.e. the space 
that remains relevant for tasks state manipulation tasks 
where we start in the state |1) when we remove all the 
CSOs from the single excitation subspace. Now we wish 
to identify the ASOs on T-La, as these are the only ones 
that could affect our state manipulation task, potentially 
reducing the achievable fidelity below that of Eqn. ©. 

Theorem 1. In a pendant- controlled system with the 
pendant quhit indexed as qubit 1, if Gq is a graph which 
is connected on vertices 2 to N with positive weights, 
d[^^ > 0, a necessary and sufficient condition for the exis- 
tence of an ASO that applies to T-La is that Gq is bipartite 
(two-colorable). 

Proof. We start by noting that if G is bipartite, we can 
easily give an ASO within the single excitation subspace. 
For a graph to be bipartite, it means that we can parti- 
tion the vertices into two sets {V}i,yB} such that edges 
only connect between the sets, not within a set. If we 
apply a product of Z operations on one side of that bi- 
partition, then every {XX -\- YY) coupling is affected 
by exactly one Z operator, with which it anticommutes. 
Hence, this is an ASO, and applies in all excitation sub- 
spaces. 

In case there are additional CSOs that apply within 
the single excitation subspace, we need to justify that 
this ASO will also apply to just Ha- Let M be the matrix 
representing the ASO in the single excitation subspace. 
We can write it as the diagonal matrix 

M=Y^ \n){n\- I")H- 

neVA neVg 

As with all ASOs, M imposes that the eigenvectors of A, 
\Xn) either have eigenvalue or arise in ±A„ pairs with 
the eigenvector of eigenvalue — A„ being given by M|A„). 
This is proven by using the anticommutation, 

A{M\Xn) - -MA\Xn) - -A„(M|A„)). 

Now let's imagine there's another CSO that applies to 
both A and G. Thus, there must be a unitary matrix that 
acts on A and G to create the subspace structure. This 
additional structure is composed of eigenvectors |A„) of 
A with (2|A„) — 0. If these eigenvectors are subject 
to the ASO M, then the remaining subspace must also 
be subject to it. If (2|A„) = 0, then is also holds that 
(2|M|A„) = 0, so for any A„ ^ 1, eigenvectors can be 
removed in pairs. These pairs manifest within A or C as 
a subspace 

5= |A„)(A„|-M|A„)(A„|Mt. 
Since = 1, we have that 

{S, M} = 0. 

On the other hand, if A„ — 0, then we just have S — 0, 
and thus {S,M} = So, for each subspace that we re- 
move, the matrix M splits into two components, an ASO 



for the removed subspace, and an ASO for the remaining 
one. 

If M were to remain diagonal under the action of that 
unitary, then M would also have the same subspace struc- 
ture, and we could just extract that "Ha component. If 
we are removing an eigenvector |A„) from the space, then 
it is because it has (2|A„) = 0, which will also be true for 
the eigenvector M|A„) i.e. for non-zero eigenvalues, the 
eigenvectors must be removed in pairs. 

Moreover, it is the only relevant ASO in the single 
excitation subspace. 

Lemma 1. Having taken into account all CSOs, T-La for 
a pendant- controlled system has no more than one ASO. 

Proof. We can take the restriction of C to Ha to be C = 
|1)(2| -h |2)(1|, and the restriction of A to Ua is A', which 
has Na eigenvectors |A„). We can always take that |Ai) = 
which is a value eigenvector of A' . 
Assume there exists an M which anticommutes with 
both G and A' . Anticommutation with G means the 
structure of M must be such that |1) and |2) are eigen- 
vectors with eigenvalues of equal magnitude but opposite 
sign. We can fix the scale factor of the matrix by specify- 
ing that 1 2) is the eigenvalue. Anticommutation with 
A' implies that eigenvectors come in ±A„ pairs, with a 
single eigenvector if Na is even. We order the eigenvec- 
tors such that the eigenvalues A„ = — A„+„^ for n = 2 to 
l + Ua with Ua — [(A^a " l)/2j . Heuce, we can write M 
as 

M = -f /3„(|A„)(-A„| + I - A„)(A„|). (11) 

n=2 

If Na is even then A' has a second 0-eigenvalue (IAat^)) 
which adds a term /3Ar„ IAa?^) (Aat^ | to M. Let U = 
be the unitary operator that diagonalizes 
A'. Then UA'W = 

f diag(0,A2,A3...Ai+„^,-A2,-A3...) Na odd 
\ diag(0, A2, A3 . . . Ai+„^, -A2, - A3 . . . , 0) Na even. 

The action of U also simplifies M, UMW = 

r (-l)©(A®diag(/32,/33..-/3l+nJ) ^aOdd 

\ (-l)©(A(g)diag(/32,/33---/3i-Hna))®/3A'a A^a even. 
As |2) is a +1 eigenvector, we have M\2) — |2). Defining 

I , r {1(B{H<»1))W\2) Na odd 
'"'^ \ {1 (B (h (8) 1)(S I) W \2) Na even 

where H is the Hadamard gate, it follows that \v) is a 
-1-1 eigenvector of {H (g) 1)M{H (8) 1), which is a diagonal 
matrix. Evidently, one of {n\v) and (n -I- na\v) must be 
for n > 2. However, both cannot be because that 
would imply {n\v) ± {n ~\- na\v) =0, and thus {n\W\2) — 
(A„|2) — 0, which is not allowed since there are no CSOs 
on "Ha. This uniquely determines M. If Na is even, the 
remaining coefficient /J^r^ = 1 as we must have {Na\v) = 
(Awj2)7^0. ° □ 
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It remains to prove that non-bipartite graphs do not 
have an ASO on Ha- Let us assume the contrary - there 
exists an ASO, M. This means that ah the eigenvec- 
tors |A„) with eigenvalue A„ 7^ have a counterpart 
M|A„) — I — A„) of eigenvalue — A„. Moreover, as argued 
above, |2) is an eigenvector of the symmetry operator, so 
these eigenvectors have the same value of a„ := (2|A„). 
Following a similar route to ^ , this assumption imposes 
that 



(2|A2'=+i|2)=^KpAf+i=0. 



(12) 



We can also expand the (2|A^'^+'^ |2) in terms of the cou- 
pling strengths dj-^"* ; it is a sum of paths of 2fc -I- 1 steps 
starting and ending on vertex 2, where for each path we 
just take the product of the coupling strengths along that 
path. The existence of such a path implies the existence 
of an odd cycle in the graph i.e. for bipartite graphs, 
^2|^2fe+i|2^ = 0. For a non-bipartite graph, however, 

there do exist odd loops. If all the df^^ > 0, then the 
product around any given path is positive, and the sum 
over all paths must also be positive i.e. there must exist a 
k for which (2| |2) > 0, and we must see such a path 

if the graph is connected and all the coupling strengths 
are positive; Eqn. p2p must be violated, and there can- 
not be an ASO [34]. Whether we operate in the single 
excitation subspace or Ha does not matter to this argu- 
ment as the only difference is the presence or absence of 
eigenvectors with a„ = 0, which therefore do not affect 
(2|A2'=+i|2). □ 

Theorem 2 (Maximum transfer fidelity). Let A corre- 
spond to the single excitation subspace of an XX Hamil- 
tonian with positive weights, and a coupling geometry 
specified by a graph Gq that is connected and not bipar- 
tite. In conjunction with a pendant-type control C , the 
maximum fidelity for transfer of an initial state |1) to a 
target state |^) in the single- excitation subspace is 

= M^Ki)\^ (13) 

{(fe,£)|(l|C|Afc.,)=0} 

and this bound is attainable, provided {\Xk.e)} is an eigen- 
basis of A with I — 1, . . . , ranfc(nfe) satisfying 

{\\C\\kx) = 0, Vfc, > 1, and (14a) 
n,„|Afe^f) =0, Vm ^ A: (14b) 

where is the projector onto the k-th eigenspace. 
Proof. Let us assume that A has an eigendeconiposition 

^-^A„|A„)(A„|. 

n 

We can always take |Ai) — |1) and Ai = as |1) is 
an eigenstate of A with eigenvalue 0. Moreover, by the 
definition of the pendant-type control system the only 



connection of this subsystem to the rest of the network 
is through state |2) via C. Define the C-overlaps 



(1|C|A„) - (2|A„), 



(15) 



and let Iq = {n|a„ ^ 0}. If a„ = then the subspace 
spanned by |A„) is decoupled and cannot be accessed. 
Hence, any state produced from the initial state using 
the control must be a superposition of eigenstates |A„) 
that have non-zero overlap with |2), 

nS/o 

Writing the target state 1-0) in the ^-eigenbasis \tp) = 
J2n (^"'^ I'^n) , the overlap of the output state with is 



As we must have \ln\^ = 1 and 7„ = for n ^ /q, 
the overlap is maximised if we choose 



(V')|2 



Vn e Iq. 



in which case we obtain 



Fl^lh — 



W|2 



nelo 



= i-Ei/3i^^i' = i-EK^i^")i'- 



A non-bipartite graph has no ASOs, so we have 
frequency-selective control of each transition |1) 1-^ |A„) 
for which the transition probability |(Ai|C|A,i)p does not 
vanish by using a Rabi oscillation of the control field C 
with frequency wi„ = A„ — Ai = A„. There are some 
special cases which we have to handle. It can be the 
case that there are two eigenvectors with eigenvalue ±A. 
Provided that the |a„| are distinct, a Rabi oscillation 
still suffices via judicious choice of pulse timing. If the 
|a„| are also equal, one has to use a different technique - 
we find a different level \k) which we populate first, and 
then perform a two frequency (Raman) transition with 
frequencies A — A and k — A, where A is a detuning. The 
only way that this technique can fail is if the system is so 
structured that all eigenvectors satisfy this pairing prop- 
erty, or have eigenvalue i.e. there is an ASO present, 
but we know there is no such ASO. 

We can therefore, in principle, transform the state |1) 
to the output state by applying a sequence of pulses 
resonant with the transition frequencies A„ as in [5|, for 
instance. With arbitrarily slow and weak pulses, the off- 
resonant excitations are arbitrarily weak, however off- 
resonant excitations could be suppressed, transfer times 



reduced, and other constraints introduced, by using op- 
timal control pulse design [l^ . In the case of an on-ofF 
control switch, we could try to match the frequencies in 
the Fourier decomposition of the square pulses with the 
resonant frequencies as pointed out in Q. □ 

In the case of bipartite graphs, there is an ASO present, 
and attempting a Rabi oscillation using C to address the 
transition between |1) and |A„) equally couples to the 
transition between |1) and | — A„), so we can only prepare 
paired eigenvectors in states 

e*^"*|A„)+e-*^"*|-A„) 

i.e. they have the same mod-square amplitude. In a sim- 
ilar way to Example 1, we will prove that this happens 
however we attempt to control the system, and give a 
physical interpretation of the states satisfying this prop- 
erty. 

Lemma 2. A state transformation task on a bipartite 
pendant controlled graph can be implemented up to the 
bounds specified in Theorem if and only if the target 
state has (up to a global phase) real amplitudes on one 
of the bipartitions, Va and imaginary amplitudes on the 
other, Vb- 

Proof. The states that we can make are of the form 
IV'targ) = ^/?„(e'^^"|A„) -f - A„)) 

n 

with real The bipartite graph in the subspace Ha 
(any weight of the target state outside this space is ex- 
plicitly forbidden by Theorem[5]) has exactly one ASO. It 
suffices for us to describe its effect on the single excitation 
subspace since we have already seen how this symmetry 
splits under the action of any further CSOs. 

Since all eigenvectors I A„) of AhaveM|A„) = |— A„), |A„) 
and I — A„) have real amplitudes 7„m with equal modulus 
on every spin to, and the relative phase in amplitude, ±1, 
is determined entirely by which side of the bipartition 
they sit on. Hence, 

)|to) 

+2i E y^/?«7nmSin(0n)|TO). 

All the amplitudes on Va are real, all those on Vb are 
imaginary and all such amplitudes are reachable. 

To see that the transformation is otherwise impossible, 
consider the operator 

Ub= U \fz'n. 
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FIG. 1: Controlled interaction between spins 1 and 2 marked 
double in gray 



If we transform Hq and He according to Ub , all {XX + 
YY) terms are transformed into {XY — YX) terms, 
which are imaginary. Hence, all evolutions —iUB{HQ + 
f{t)Hc)Ug are real, and the corresponding unitary op- 
erators are real. By starting with a state J7b|1) = 
which is real, the amplitudes of all output states Ub\'4>) 
must be real. Thus {ip) must be real on spins that are in 
Va, and imaginary on those in Vg. □ 

In particular, this means that state transfer tasks are 
unaffected by ASOs. The phase result makes a lot of 
sense because we know that evolution under a Hamilto- 
nian i {XX + YY) leads to a hopping term between two 
vertices, but introduces a phase factor of —i for each hop. 
Hence, on a bipartite graph, if we start on a single site 
with a real amplitude, all those that are an even number 
of hops away (which is only well defined on a bipartite 
graph) must have a real amplitude, whereas those that 
are an odd number of hops away must have an imaginary 
amplitude. 

Alternatively, the limitation is easily overcome - we 
can introduce a second control field C = «|1)(2| — i|2)(l|, 
which corresponds to a coupling ^{XiY2 — F1X2). As 
this coupling is not real, it does not satisfy ([3]), although 
it does still anticommute with the symmetry operator. 
Using a Rabi oscillation of cos(aji„i)C ± sin(cji„t)C" al- 
lows us to selectively make the transition of one or other 
of the eigenvectors. This illustrates the important dif- 
ference between the symmetry condition of Eqn. ([3]) and 
anticommutation if the Hamiltonian is not real-valued. 
With this additional control in place, the conditions of 
Thm. [2] are all achievable. However, we will not further 
consider this control. 

In summary, we see that in all cases the maximum pos- 
sible state transfer fidelity is entirely determined by the 
energy eigenstates of A which have zero overlap with |2), 
the point where the pendant control meets the rest of the 
system. Entanglement generation, however, is addition- 
ally affected in bipartite graphs by whether the pair of 
qubits that are to be entangled are separated by an even 
or odd number of edges, and the relative phase to be gen- 
erated between the pair - we can either make \n) + \m) 
for even distances, or \n) + i\m) for odd distances. 
Example 2: First we examine a uniformly XX-coupled 
iV-spin chain with a fork at the end (see Fig. [T]), initial- 
ized with a single excitation localized on the first qubit. 
Diagonalizing the adjacency matrix shows that the only 
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FIG. 2: Controlled interaction between spins 1 and 2 marked 
double in gray 

eigenstate with zero C-overlap is |A_) ~ ~ I '''))■ 

Consider state transfer to |6); when we write this in the 
adjacency matrix eigenbasis, the coefficient of the |A_) 
eigenstate is Z?!*'^ = l/'s/Z, giving a maximum transfer 

fidelity Fi^g = 1 - l/^i^^ = 1 - \^\^ = \\ the same 
result is obtained for transfer to |7) by symmetry. This 
is to be expected; the two end qubits are in some sense 
indistinguishable, and we cannot direct information flow 
to either one specifically. At all times, the control ad- 
dresses both simultaneously, and so we obtain a super- 
position state over these two qubits. This is made clear 
in the algorithmic computations. In this case, the Hamil- 
tonian is symmetric under the change of basis CSO J^-j 
taking J67|6) = 1 7), J67|7) = |6), Jgyln) = |n) other- 
wise. Performing the computations above shows that the 
Hamiltonian in this subspace generates a dynamical Lie 
algebra of so (6) ® 1, where the 1-term corresponds to 
the subspace spanned by the dark state -^(|6) — |7)). If 
our initial state has no overlap with this state then it 
remains unpopulated and "dark" throughout. Alterna- 
tively if our initial state is contained in this subspace, it 
will be preserved throughout the system's evolution. On 
the other hand, we can choose to produce states such as 
( 1 6) -I- 1 7) ) / -\/2, which is an entangled state, or states such 
as 1 5), which would allow perfect transfer of an unknown 
qubit state from qubit 1 to qubit 5. However, without 
the additional control C", it is impossible to perfectly 
produce states such as "^(|3) -|- |4)) as a result of the 
bipartite nature of the graph, and the disruption caused 
by the ASO. 

The results in the previous example are intuitive, but, 
for other systems, identifying the symmetries and dark 
subspaces may be less intuitive. 

Example 3: If we make a small modification to the 
system in Fig. [T] by moving the pendant vertex as shown 
in Fig. [21 then additional symmetries and dark states 
arise. Since the graph is bipartite, with partitions V\ — 
{1, 3, 4, 6, 7} and V2 — {2, 5}, we can immediately observe 
that there is an ASO 

Af = l-2(|2)(2| + |5)(5|), 

derived from .^2-^5, which anticommutes with the Hamil- 
tonian. We again observe the permutation CSO 

J67 = l-(|6)-|7»((6|-(7|), 



which commutes with the Hamiltonian. These two sym- 
metry operators give rise to an eigenspace of eigenvalue 
zero which has no overlap with 1 2) and is spanned by the 
vectors 

|Ao,i) = ^(|6> - |7)), 

|Ao,2> = ^(2|3)-2|4) + |6) + |7)). 

No other eigenvectors have zero overlap with |2) thus the 
maximum transfer fidelity to |3), for example, is F\y^^ = 

1 - - = f as /3f) = and fif = ^. 

Similarly, the state -^(|6) -|- |7)) can only be produced 

with fidelity 1 - i = |. 

IV. CATALYTIC EXCITATIONS AND THE 
CLASSIFICATION OF DARK STATES 

So far we have discussed access to states in the first- 
excitation subspace using paths entirely within this sub- 
space, as assumed universally in prior works. The con- 
trolling interaction was assumed to commute with the 
total spin operator Sz and was thus unable to produce 
paths which move between excitation subspaces. To ac- 
cess dark states, we could modify the control so that 
it does not commute with Sz^ breaking the symmetry 
and thus the excitation subspace structure. Doing so 
could lead to difficulties: the resulting subspaces will be 
of much larger dimension, and this will increase the com- 
putational difficulty of finding controls. The higher di- 
mensional systems are more likely to contain uncharac- 
terized resonances which makes it harder to produce con- 
trols which confine the excitations to the desired path in 
state space. However, we shall now show that moving to 
higher excitation subspaces can sometimes be advanta- 
geous and need not complicate the control too much. 

In general, breaking symmetries requires changing the 
control Hamiltonians or adding new controls. However, 
in some cases we can avoid certain symmetries by moving 
to higher excitation subspaces, using our ability to intro- 
duce additional excitations on the first spin, which has 
been implicitly present as a control within the initializa- 
tion stage, but we have not explicitly considered it when 
deriving accessible states. To produce a desired target 
state l/fe) in the fc-excitation subspace, we may be able 
to bypass symmetry-induced restrictions in this subspace 
by temporarily introducing additional excitations, which 
we shall call catalytic excitations. To do so, we pass to 
the Z-excitation subspace which contains a state of the 
form l/fe) (g) |(/)) in its reachable set, and finally extract 
the excess excitations to leave Not all symmetries 

can be avoided this way; permutation symmetries which 
do not act on qubit 2 certainly cannot be avoided. How- 
ever, in many cases, adding a single extra excitation will 
suffice. 

From [i^l, we have 
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Definition. The matrix representing the second excita- 
tion subspace of Hq is given by 

i |((fc, Z| - (/, fc|) [In^A + A^ In] j) - \j, (16) 

for m := {N — l)k + I and n := {N — l)i + j, where \i,j) 
denotes the vertices at which the excitations are localized. 

Similarly, (C) is the action of the control upon this 
subspace. The connectedness of this graph under ba- 
sis changes (and thus the irreducibility of the dynamical 
Lie algebra) is not easily determined from the connect- 
edness of A and C [s^. In other words, the existence, or 
lack thereof, of ASOs and CSOs other than permutation 
symmetries in one excitation subspace does not necessar- 
ily imply their existence in another excitation subspace, 
therefore we must search for them in each subspace sepa- 
rately. It is this feature that potentially enables catalysis, 
illustrating that, in general, it is not sufficient for anal- 
yses of state transfer in spin networks to restrict to the 
single excitation subspace. 

Theorem 3 (Sufhcient condition for increased control 
through catalytic excitations). A system has an accessi- 
ble pair (IV'm), IV'out)) 'in the single excitation subspace 
using catalytic excitations (assuming {l\^pout) — ^) if 
it has accessible pairs (IV'm)? IV'')) '^^d (| 1, '0'), |1, V'out)) 
with (IIV'') = 0. 

Proof. Initialize the system in the state |1) and trans- 
form it to the state |^') in the single-excitation subspace, 
which by hypothesis is accessible and satisfies (llV'') = 0. 
Introducing a second catalytic excitation at the first site 
then produces the state |1) ® \tp') in the second excita- 
tion subspace, and by hypothesis, we can create the state 
|l)i8)|V'out)- Removing the catalytic excitation on the first 
site, we are left with our desired target state, IV'out)- CH 

We now give an example to prove that such a situation 
can arise, even when (|i/'in), IV'out)) are not an accessible 
pair when confined to evolution solely within the single 
excitation subspace. 

Example 4: For the system given in Fig. [2l we saw in 
Example 3 that the target state |3) is not perfectly acces- 
sible from the initial state |1). As before, we look at the 
eigenbasis of A where we have eigenvectors correspond- 



1,3 ^1,2 \A \,5 ^1,+ 



ing to eigenvalues of Aq = and A±^± — ±y (5 ± V5)/2. 
Labelling the corresponding eigenvectors |Ao) and IA4-+), 
I A |_), etc. we have 



(|A++) + |A_+)) 
(|A+_) + |A__)) 



|3> = ^/5lo) 



10- V2 

T0 + V2 




O 



viy 3,- 2,- 4,- 5,- 



FIG. 3: Controlled interactions marked double in gray; i,j 
represents product state \i)\j); 4- represents -^(16) + |7)), - 
represents -^(|6) — |7)). Note that the controlled couplings 
are simultaneously controlled, not independently. 



and starting from a|0) -I- /3|1), the state with maximum 
overlap with |3) we can achieve is 



|3)- 



|4> + Viij(|6) + |7))) 



with maximum fidelity [36(1 



Instead of this, we could first produce 
^+V2 



+ 



4^/5 



(|A+-> + |A-)) 
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(7-V5)(|A++) + |A_+)) 



defining j?/;'), and then introduce a single catalyst excita- 
tion at spin 1, 



all) 



10 + V2, 



4\/5 



l,A+_) + |l,A__)) 



4^/5 



+ \/y^(7-V5)(|1,A++) + |1,A_+)) 



We have moved into a superposition of the 1 and 2- 
excitation subspaces, and our control will keep the state 
in these subspaces. To show that we can reach the target 
state using catalytic excitations we need to prove that 
the set of states 

{(|l,V'),|l,V'out)),(|l),|l))} 

are simultaneously accessible such that the |0) compo- 
nent of the unknown state is also correctly transferred. 
This 2-excitation subspace has no dark states which are 
left uncoupled to the controlling interaction, except those 
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which are a product of ;^(|6) — |7)) and a single exci- 
tation on one of qubits 1 to 5. Examining the graph 
produced by the Hamiltonian for the 2-excitation sub- 
space (Fig. [31 formed as in [13]), we can see that there 
are no eigenstates of the Hamiltonian for the left-hand 
component (i.e. the connected component to our initial 
state) that have no overlap with all the states our con- 
trol connects to. This is the analogue to checking which 
states had (1|C|A„) — above. 

We select an eigenvector |A) in the second excitation 
subspace. Based on this choice, we apply a field of 
di2 oc cos((A++ — A)i) on the controlled coupling between 
spins 1 and 2, and, neglecting off-resonant excitations, we 
obtain an effective coupling 

77icx|l,A++)(A| + |A)(l,A++| 

+ |l,A_+)(-A| + |-A)(l,A_+| 

which can transform the state to 



+ I-A)) + ^^1^-^(11, A++) + |1, A_+)) 



Similarly, neglecting off-resonant excitations, applying a 
field Ji2 (X cos(Ai) on the controlled transition gives rise 
to an effective Hamiltonian 

iJ2 (X |1, Ao)(A| + |A)(1, Aol + |1, Ao)(-A| + | - A)(l, Ao| 
which can transform the previous state into 



^^(|l,A,_) + |l,A--)) + /i(|A) 



4^/5 



-(|l,A+_> + |l,A__» 



+ y||l, Ao> + "^^"^ (11, A++) + |1, A_+» 



«|1) 



-/3|1,3). 



state ^ 



Removing the catalyst excitation from the first spin gives 
the desired final state, a|0) + /?|3). Introducing this ex- 
citation allowed us to move between different subspaces 
induced by the symmetry. 

Notice that the catalyst excitation cannot make the 
— 1 7)) accessible from an initial state |1) be- 
cause the graph derived from the Hamiltonian will have 
two components, with one component comprising states 
of the form !■(/;) ® — 1 7)) in all the excitation sub- 

spaces. Therefore, we can say that there are two types 
of dark state, one of which one is not as dark as the 
other. Dark states arising from a permutation of the 
qubits [17[ are truly dark - they induce a splitting of the 
Hamiltonian which persists in all excitation subspaces, 
and which cannot be overcome without introducing a 
symmetry-breaking control. The second class of dark 
states are artifacts of the particular eigenvectors found 
in each subgraph, whose underlying symmetries may be 



bypassed in higher excitation subspaces, potentially ren- 
dering them accessible. However, it should be noted that 
if A is a simple one-dimensional chain then these symme- 
tries do persist, as can be proven via the Jordan- Wigner 
transformation [5|]. 



V. SYSTEM IDENTIFICATION 

Our control method in Sec. IIIIl based on geometric 
control, is sub-optimal and, if the system Hamiltonian is 
known, there are various ways to derive control pulses 
that saturate the bounds. However, these require some 
basic knowledge of the Hamiltonian. In particular, given 
the promise of pendant control, we only need to know 
the eigenvalues A„, «„ (the overlaps of the eigenvector 
with 1 2), the single excitation on spin 2), and some in- 
formation about the target state ~ its decomposition in 
terms of the eigenvectors. This decomposition must be 
specified with a particular phase convention for the eigen- 
vectors, which we shall choose such that > for all 
n. Numerical simulations further show that it is possible 
to derive pulses to achieve this threshold even when the 
Hamiltonian of the system has not been determined, us- 
ing optimization over the record of multiple single-basis 
projective measurements on the target qubit, as devel- 
oped in [9)] . Indeed, we shall now prove that it is possible 
to identify the control parameters A„ and a„ so that, if 
given a target state expressed in terms of the eigenvec- 
tor decomposition, the state transformation task can be 
achieved without any further promises on the system. 

For this, consider the same drift Hamiltonian A and 
control field C that we have used so far within the single 
excitation subspace. By determining the eigenvalues An 
and overlaps a„ = (2|A„) of A', the restriction of A to 
the space "Ha, then if we are given a state transformation 
task starting from |1), with an output state described in 
terms of the amplitudes of different eigenvectors, we can 
implement that transformation optimally. Of course, in 
order to translate from a specification in terms of the spin 
basis will require further promises about the structure of 
A. To date [2l|,[23|, such assumptions represent a massive 
restriction, so it is interesting to identify what aspect of 
the protocol it is that requires such assumptions. 

Let us define 

i? = ^ + eC, 

and assume that H has eigenvectors \r\n). Our level of 
control can certainly allow us to measure, as a function 
of time, quantities such as 



(l|e- 



-im 



11)1^ = 



a Fourier transform of which will reveal the eigenvalue 
differences Irjn — rjml and overlaps |(l|77„)p. By assuming 
e to be small, we can perform a perturbative expansion. 
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expressing these quantities in terms of the properties of 
A' that we wish to determine. In particular, for n 1 
and for A„ 7^ 0, we can write that, to first order, ry„ = A„ 
and 



I '7" 



= IA„ 



g(l|C|A„) 



ID 



i.e., |(l|77n)| = ean/Xn, so the a„ and |A„ - A^l are 
readily extracted for any eigenvector which does not have 
a„ = i.e. as one would expect, we cannot learn the 
properties of the states outside the space Ha- If there 
is an eigenvector with Am — and am 7^ 0, then we 
have to apply degenerate perturbation theory, but this 
just means the two energy levels (including |1)) mix, 
with energies ±eam- Again, the relevant quantities can 
be extracted. 

From the differences |A„ — Am|, the eigenvalues can 
only be identified up to an overall global shift, Ai, and 
the difference between A„ — Ai and — (A„ — Ai) cannot 
be distinguished. The global shift only alters the global 
phase in the subsequent state transformation protocol, 
and is therefore largely irrelevant, but we can still ar- 
gue how to determine both of these features. Firstly, 
observe that neither affect a system with an ASO - we 
know that H = —H. Hence, we only have to deal with 
the case where the levels can, in principle, be addressed. 
For a given eigenvector |A„), we know q;„, but not the 
eigenvalue. By determining two such values, all ambigu- 
ity can be removed. If we have an eigenvector |A„) with 
eigenvalue A„, then provided there is not also an eigen- 
vector with —X„ (if there is, we can adapt for it), we 
can perform a Rabi oscillation of C with frequency A„, 
and in some fixed time, as specified by |a„p, we return 
to our initial state, having gone down to a probability 
of |a„p of being found in the state |1) half way through 
the protocol. Hence, for a unique value of janP, we can 
scan through all possible frequencies and find the relevant 
value of I A„ I . In order to discover the sign of this energy 
(since a Rabi oscillation addresses ±A„ simultaneously), 
we instead prepare |0) + e"^|l) for various values of <j>, 
and transform it to |0) + e^^jAn), leave it for some time t 
so that it evolves to |0) -I- e'*^''^"'"^"*'' |A„), and then return 
it to the initial state |0) -I- e^^'^~^^"*^l) . By measuring 
the probability p = cos^(|((/) -I- A„i)) of getting the |-f) 
measurement result, then provided |A„|t < tt, the sign of 



dp 



6=0 



is the same as that of A„ . 

This analysis reveals that within a pendant controlled 
system, it is, in principle, possible to extract the values 
A„ and that we need to know, although in practice 
one would undoubtedly develop more sophisticated pro- 
tocols which do not rely on measuring extremely small 
quantities. Consequently, we observe that to implement 
a complete identification protocol such that we can pro- 
duce a state described in the position basis rather than 



the eigenbasis of A, it is necessary to be given sufficient 
information on the system to be able to derive the eigen- 
vectors from the values An and an ■ A chain is a particu- 
larly natural candidate and can be extended to chain-like 
systems [21j, 22^ , although one could come up with many 
other variants 



VI. CONCLUSIONS 

In this paper, we have demonstrated how the capabil- 
ity of a spin network for state transfer can be computed, 
which has obvious potential bearing on the design of 
quantum information routers. All CSOs produce states 
which are robust to dephasing noise on the pendant con- 
trol spins, and those for which catalysis fails (such as per- 
mutations) produce states which are robust to all noise 
on the pendant control spins. The difference between 
the two can be interpreted in terms of selection rules 
[2^ : the selection rule preventing population of some 
states can be overcome by shifting to a different excita- 
tion subspace, but superselection rules prevent access to 
the truly dark states. This could direct the development 
of spin network memories in terms of decoherence-free 
subspaces, i.e. the dark states produced by the permuta- 
tion symmetries of the excitation subspace graph could 
be used to store states. Truly dark states will be robust 
to the introduction of extra excitations in the system, 
but will need a corresponding permutation symmetry- 
breaking local control to enable their initialization and 
collection. The weaker class of dark states will be im- 
mediately vulnerable to the introduction of excitations 
elsewhere in the system; nevertheless they will be bet- 
ter protected than non-dark states and correspondingly 
easier to access using local controls than the truly dark 
states. The suitability of each method of storage will 
depend on the architecture employed and the relative 
importance placed on access and permanence. This is 
somewhat analogous to the relative advantages of RAM 
and ROM in conventional silicon computing. 

We have seen how, in the pendant controlled spin net- 
works, the necessary and sufficient condition for the ex- 
istence of an ASO (on T-la) is that the underlying graph 
should be bipartite. This link with the underlying graph 
structure is extremely interesting because it is much less 
artificial than previous studies which have found such a 
relation [13] . These studies imposed a direct relation be- 
tween the Hamiltonian and the adjacency matrix of the 
graph by setting them equal, whereas our only constraint 
is that the graph structure specifies where coupling coef- 
ficients take on a non-zero value. 

In the future it will be interesting to see when the 
state with optimal overlap can be achieved for arbitrary 
Hamiltonians, not just those with pendant controllers. 
A trivial extension of the present study is one where the 
graph A is composed of two connected components, and 
the control C couples between a single vertex of each of 
the two components. Further consequences of the control 
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required to introduce a state to transfer should also be 
examined [30]. 
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Appendix A: Identification of Symmetries 

Given the Hamiltonian Hq of an iV-dimensional quan- 
tum system and its M control Hamiltonians i/m, it is 
relatively easy to determine the symmetries it possesses 
and thus limit the possible operations that can be real- 
ized. Eirst we check for decomposability into a block- 
diagonal structure (i.e. we find the CSOs). This can be 
done by transforming each of the M Hamiltonians Hm 
from its iV X matrix representation into its Liouville 
representation, where it is represented by an N'^ x A^ 
matrix i^™) with entries 

^s!r^ '■= {Hm)nM^j,e ^ {Hm)j,tSn,k (Al) 

for r ^ [j - l)N + k and s ^ {£ - 1)N + n. Repeat- 
ing this for all system and control Hamiltonians and 



stacking the resulting matrices L*^™^ vertically gives an 
(M -I- 1)A^ X A^ matrix. Any standard routine can be 
used to calculate the null space of this matrix, and the 
column vectors of length A^ spanning the null space, 
"unstacked" into a A x A matrix, to give the CSOs J 
that simultaneously commute with all -ffm, [H^m J] = 0. 
This is a conceptually very simple approach. More effi- 
cient alternatives for simultaneous block-diagonalization 
such as the algorithm in p4| . which requires handling 
matrices of lower dimension, exist. 

To find the ASOs, we go through much the same pro- 
cedure: for each indecomposable subspace, indexed by d, 
of each Hamiltonian Hm, we take the complexified sys- 
tem and control Hamiltonians restricted to the subspace, 
iHm^d, remove the trace 

Hm.d iHm,d - Tr{iH,n,d)-^, (A2) 
and compute the Liouvillian representation 

L^™''') := {Hm,d)n,kSjJ + {Hm,d) j,l5n,k- (A3) 

As before, we do this for all system and control Hamil- 
tonians, stack the resulting Liouville operators, and cal- 
culate the null space of the resulting matrix. Unstack- 
ing the null vectors again results in ASOs J such that 
Hm,dJ + JHm,d = for aU Hm,d. The ASOs J define the 
possible superpositions within the subspaces spanned by 
the eigenvectors of the J CSOs. The ei genv alues of the 
J may help identify the Lie subalgebras (25l [26j . 
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